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Abstract. Using spectra of normal emission line galaxies from the First Data Release of the Sloan Digital Sky Survey (SDSS) 
we have investigated the relations between the extinction C(H/3) as derived from the Ha/FJ/3 emission line ratio and various 
global parameters of the galaxies. Our main findings are that: 1) C(H/3) is linked with the galaxy spectral type and colour, 
decreasing from early- to late-type spirals. 2) C(FJ/3) increases with increasing metallicity. 3) C(H/3) is larger in galaxies with an 
older stellar population. 4) C(H/3) is larger for more luminous galaxies. 5) The extinction of the stellar light is correlated with 
both the extinction of the nebular light and the intrinsic galaxy colours. We propose phenomenological interpretations of our 
empirical results. We have also cross-correlated our sample of SDSS galaxies with the IRAS data base. Due to the lower redshift 
limit imposed to our sample and to the detection limit of IRAS, such a procedure selected only luminous infrared galaxies. We 
found that correlations that were shown by other authors to occur between optical and infrared properties of galaxies disappear 
when restricting the sample to luminous infrared galaxies. We also found that the optical properties of the luminous infrared 
galaxies in our SDSS sample are very similar to those of our entire sample of SDSS galaxies. This may be explained by the 
IRAS luminosity of the galaxies originating in the regions that formed massive stars less than 1 Myr ago, while the opacity 
of galaxies as derived from the HafHf3 emission line ratio is due to diffuse dust. We show some implications of our empirical 
results on the determination of global star formation rates and total stellar masses in normal galaxies. 

Key words. Galaxies: spirals - Galaxies: abundances - Galaxies: evolution - Galaxies: ISM - Galaxies: stellar content - 
Galaxies: ISM: dust, extinction 
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Bianchi et al. JT9961 . and Witt & Gordon (HQ). Infrared data 
from the IRAS survey provide detection of dust through its 
emission in the infrared. However, the interpretation of the in- 
frared dust emission in terms of dust content is not straightfor- 
ward (see e.g. Sauvage & Thuan lT^ Dale & Helou l2UU2l . A 
different approach is to study the transparency of galaxies with 
respect to background light (see Keel & White I200T1 and refer- 
ences therein for methods using one nearby background galaxy, 
and Gonzalez et al. 1998 for methods using distant galaxies as 
background sources). Unfortunately, this method suffers from 
extremely low statistics. 

Another extinction indicator that can be used in galaxies 
presenting emission lines is the observed Balmer decrement. 
As is known, the intrinsic intensity ratios of hydrogen recom- 
bination lines in ionized nebulae have a negligible dependence 
on the emission conditions (Osterbrock 1989 1 so that the ob- 
served values are a consequence of the dust extinction be- 
ing different at the relevant wavelengths. This has for exam- 
ple been used in the influential work by Calzetti et al. {1994) 
to determine an empirical extinction law for starburst galax- 
ies. However, few studies have used this method to investi- 
gate the relation between extinction and overall galaxy type 
(e.g. Hubble type). In order to do this, one needs spectra of 
galaxies spanning the entire range of Hubble types. The atlas 
of Kennicutt ( 1992 > provided such a data base for a limited 
sample of galaxies. Among these galaxies, 15 were considered 
by Sodre & Stasinska ( 1999 hereinafter SS99) to be normal 
emission line galaxies, and these authors showed that the ex- 
tinction at HB as measured by the HafHB ratio (in the remain- 
ing of the paper we will call it the Balmer extinction) decreases 
steadily from early-type to late-type spirals. Later, using the 
Nearby Field Galaxy Survey (NFGS) of Jansen et al. (2000a 
2000b), which provided adequate data for about 100 galaxies, 
Stasinska & Sodre (2001 , hereinafter SS01) showed that redder 
galaxies have larger Balmer extinction. The Sloan Digital Sky 
Survey (SDSS), which aims at obtaining spectra of 10 6 galax- 
ies in the nearby Universe, provides a wonderful opportunity to 
study this question in more detail. 

The present paper makes use of the observations from 
the First Data Release (DR1; Abazajian et al. 2003, see also 
Stoughton et al. 2002) of the SDSS to study the relation be- 
tween the Balmer extinction and other global properties of the 
galaxies. Section 2 explains the selection of the observational 
sample and quantities used in the analysis. Section 3 presents 
the relation between Balmer extinction and other galaxy pa- 
rameters. Section 4 presents an interpretation of our results. In 
Section 5 we cross-correlate the galaxies of our sample with 
galaxies detected by IRAS, looking for additional clues on the 
origin of the observed extinction. Section 6 summarizes the 
main conclusions of this work and outlines some prospects. 

2. The observational sample 

2.1. Measurements of emission line parameters 

The equivalent widths of Ho- and HB, EW(Ha) and EW(HyS) 
respectively, measured by the SDSS team and available in 
the survey database, do not account for underlying absorp- 
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Fig. 1. Examples of two fits using our code to measure the 
equivalent widths and fluxes of HB (top panels) and Ha (bottom 
panels). The left and right figures are for spectra with signal- 
to-noise ratio in g-band of 5.9 and 12.4, respectively. We also 
show the equivalent widths (EW) in A and signal-to-noise ratio 
(S/N) of the measured emission lines. 



tion, which may be significant for objects with faint emission 
lines. Consequently, we have developed a code for measuring 
simultaneously the emission and absorption features of these 
Balmer lines (see Mateus & Sodre 2004 for additional details). 
For HB the code fits two Gaussian functions, one in emission 
and the other in absorption, allowing for different centroids and 
widths for these features. For Ha we fit, additionally, two other 
Gaussian functions to the emission lines of [N n]/i6548 and 
[N ii]/16584. The continuum is estimated through a robust lin- 
ear adjustment considering two wavelength regions on the sides 
of the lines. In Fig.[T]we show examples of the fits obtained for 
HB (top panels) and Ha (bottom panels) for two spectra with a 
signal-to-noise ratio in the g-band of 5.9 and 12.4, respectively 
from left to right of the figure. Note that these are among the 
worst cases in the sample selected for this study (see Section 
2.2). This procedure allowed us to determine, for each spec- 
trum, the fluxes and equivalent widths resulting from the ab- 
sorption and emission features of Ha and HB, as well as their 
errors and a signal-to-noise ratio (S/N) for each fit. 

We have made many simulations to verify the reliability of 
these fits and the presence of any bias in the results. We sim- 
ulated absorption plus emission lines, with realistic noise, and 
compared the input line parameters with the outcome of our fit- 
ting software. The initial parameters (essentially the equivalent 
widths of the components in emission and absorption of Ha 
and HB lines) used in these simulations were based on the re- 
sults of the spectrophotometric model of Barbara & Poggianti 



( 1997 ), for different spectral types. For each out of 5 types, we 
simulated 1000 spectral regions at Ha and HB lines, varying 
the amplitude of the noise, as well as the linear spectral disper- 
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Fig. 2. Results of the simulations for a spectral dispersion of 
1.0 A and typical noise. The input equivalent widths are based 
on the results of the spectrophotometric model of Barbara & 
Poggianti ( 1997 1 for different spectral types: Sa, Sb, Sc, Sd and 
'Extreme' case with constantly increasing star formation rates. 
These equivalent widths are compared with the values recov- 
ered by our code. The large squares are the mean values of the 
output equivalent widths, with the error bars showing their dis- 
persion. 



sion (0.5, 1 .0 and 2.0 A). Then, we applied our code to compare 
the equivalent widths of the Ha and Hfi features. In Fig. [2] we 
show a typical result obtained from our simulations, where the 
initial equivalent widths of Ha and H/3 (input parameters) are 
compared with their respective values measured by our code. 
We concluded that for S/N high enough (larger than ~ 10) the 
results recovered by our software are, within the errors, consis- 
tent with the input values and do not depend on the dispersion 
we used. 

In this study we use our own measurements of the equiv- 
alent widths of the hydrogen emission lines, and the measure- 
ments available in the SDSS database for the remaining lines. 

2.2. The selection procedure 

The spectroscopic data analyzed in this work were extracted 
from the SDSS main galaxy sample available on the First Data 
Release (Abazajian et al. 2003 1. This flux-limited sample con- 
sists of galaxies with r-band Petrosian magnitudes r < 17.77 
and r-band Petrosian half-light surface brightnesses fiso ^ 
24.5 mag arcsec~ 2 , comprising a total of 113199 galaxies. We 
first selected all galaxies with spectra obtained by the survey 
and with a redshift confidence > 0.7. This selection resulted in 
a first sample with 1 10913 galaxies. Following the conclusions 
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Fig. 3. Emission line diagnostic diagram to separate normal 
galaxies from galaxies hosting an AGN. The points represent 
the 10822 galaxies from the SDSS that have been selected ac- 
cording to the criteria described in Sec. 2.1. Solid line: the 
curve of Heckman & Kauffmann (2003); all the galaxies be- 
low that curve (9840 objects) constitute our sample of normal 
galaxies. Dot-dashed line: the theoretical curve of Kewley et al. 
(200 1) above which galaxies are dominated by an AGN; dashed 
line: empirical lower bound of normal galaxies (see Sect. 2.2). 



properties of the galaxies. We have also rejected galaxies with 
an average spectroscopic signal-to-noise ratio smaller than 5 in 
the photometric g passband. These criteria result in a selection 
of 89174 spectra of galaxies. 

For these objects, we have applied our code to measure 
the fluxes and equivalent widths of the emission Balmer lines, 
obtaining 40672 spectra with these parameters measured with 
S/N > 1 . Following the previous discussion on the Ha and H/3 
measurements, we consider only galaxies with high S/N for 
their emission lines, equal to or larger than 10, and measured 
EW(Ha) and EW(Hy6) larger than 1 A. For galaxies with strong 
emission in the Balmer lines (where the absorption equivalent 
width is much smaller than the emission equivalent width), 
our emission equivalent widths are very similar to those in 
the SDSS database. With these criteria, the sample comprises 
1 1066 spectra. In the case of duplicate spectra, we removed the 
ones with lower S/N in the g band. The resulting sample con- 
tains 10854 objects. 

We then excluded from the sample the galaxies whose 
spectra are contaminated by non-stellar activity. For this, we 
used the [O m] A5001/R/3 vs. [N n] i6584/Ha diagram intro- 
duced by Baldwin, Phillips & Terlevich (1981 1 and Veilleux 



& Osterbrock ( 19871. The resulting diagram for our sample is 
shown in Fig. [3] where we only plot the objects with measured 
fluxes for all emission lines (10822 objects). The dot-dashed 



line is a limit defined by Kewley et al. (2001 1, above which 
an active galactic nucleus (AGN) dominates the emission-line 
spectrum. The solid curve shows the empirical limit proposed 



of Zaritsky, Zabludoff & Willick ( 1995 1, we have used a red- 
shift limit of z > 0.05 to minimize the effects of the aperture 
bias in our sample, since we are interested in the global spectral 



by Heckman & Kauffmann (2003 1, below which star forma- 
tion dominates the emission-line spectrum. Galaxies below this 
curve will be referred to as normal galaxies. The dashed line de- 
limits the narrow strip where normal galaxies are found in this 
diagram. Note that this strip has roughly the same location as 
the sequence drawn by individual giant H n regions in galaxies 
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Fig. 4. Histograms of general properties of our sample of nor- 
mal galaxies. 
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Fig. 5. Some global properties of the galaxies of our sample 
as a function of their "spectral type". The top number on the 
right is the total number of data points (including a few indi- 
vidual points that may lie outside the panel frame). Below are: 
the value of Kendall's r^, then the associated probability and 
finally the Spearman rank correlation coefficient r$ (see Sect. 
2.3). 



(Mc Call et al. [T985l van Zee et al. 1998 1. Thus, in our selec- 
tion we can distinguish 9840 normal galaxies and 982 AGNs. 
In the following, our sample will contain only these 9840 nor- 
mal emission line galaxies, to avoid any interference of active 
nuclei in our interpretation of the data. 

2.3. Overall properties of our sample of normal 
galaxies 

With 9840 objects, we have over 80 times more galaxies at our 
disposal than in the SS01 study based on the NFGS sample 
(Jansen et al. 2000a 2000b I, and over 600 times more than in 



the SS99 study based on the Kennicutt ( 1992 1 atlas. Note also 
that at the spectral resolution of the SDSS, which is about 3 A, 
it is much easier to separate emission from absorption at H/3 



and Ho- than in the NFGS (spectral resolution 6 A) and in the 
Kennicutt ( 19921 atlas (spectral resolution 12-25 A). Therefore 
we are able to explore with confidence galaxies with quite weak 
H/3 emission. 

The global characteristics of the selected sample are shown 
in Fig. 4. Panel a shows the histogram of the galaxy spec- 
tral classification parameter eClass as provided by the SDSS 
for DR1. This spectral classification parameter is based on a 
principal component analysis method (Connolly et al. j 1995 
Connolly & Szalay 119991 . Note that the spectral classifica- 



tion of galaxies by Connolly et al. QT995 1 is slightly different 
from the one used in SS99 who removed from the analysis the 
wavelength regions containing emission lines. We will never- 
theless refer to the parameter eClass as to a "galaxy spectral 
type". Note that eClass is strongly correlated with the galaxy 
colour (m - z). The values of eClass, range from about -0.2 to 
0.6 from early- to late-type in our sample of galaxies. Panel 
b shows the distribution of redshifts z. Some SDSS spectra 
have residual sky lines that may affect our Ha and H/3 mea- 
surements. Despite the mask array contained in each spectrum 
file that can be used to remove any inconvenient features, like 
bright sky lines, we have chosen an alternative way by exclud- 
ing all objects in which the bright sky lines at 5578.5 A and 
7246.0 A fall on the Ha and H/3 line regions, which explains 
the bizarre aspect of the redshift distribution of this sample. 
Panel c shows the distribution of galaxy absolute magnitudes 
in the SDSS photometric r band, M(r). This magnitude is uti- 
lized because K-corrections are modest, the radiation is pro- 
duced mainly by the older stars that dominate the stellar mass, 
and uncertainties in Galactic reddening make little difference 



to the inferred galaxy magnitude (Strauss et al. 2002). To com- 
pute M(r) from the observed r magnitude given by the SDSS, 
we assumed a standard cosmology with Q.m - 0.3, = 0.7, 
and Ho = 75 km s Mpc~' and applied a K-correction (com- 
puted with kcorrect vl.ll; Blanton et al. 2003 ). 

Because of our restrictions on the quality of the spectra, we 
miss the less luminous galaxies of the SDSS (compare e.g. the 
histogram of the M(r) values in Blanton et al. 2001 1. Finally, 
panel d shows the distribution of the galaxy surface bright- 
ness, S B (in mag arcsec -2 ), defined within the radius contain- 
ing 50% of the Petrosian flux for each band, petroR50. Note 
that for this study, we do not work with a complete volume- or 
magnitude- limited sample, so that the distribution of galaxies 
among the various classes does not correspond to the true dis- 
tribution in the local universe. Also, note that, because of the 
conditions imposed on the quality of the spectra, the redshifts 
in our sample do not exceed a value of 0.25. 

Some of the overall characteristics of galaxies are in fact 
correlated. For example, galaxies of earlier spectral types tend 
to be more luminous, as shown in panel a of Fig. 5 which plots 
M{r) vs. eClass. The top number on the right corresponds 
to the total number of data points (including a few individual 
points that may lie outside the panel frame). Next, are given 
the results of non parametrical tests: from top to bottom, the 
value of Kendall's Tk, the associated probability (a large value 
of this probability indicates no correlation) and the Spearman 
rank correlation coefficient rs (see Press et al.[1992 1. Thus, we 
see that indeed M(r) is correlated with eClass. On the other 
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Fig. 6. The relation between log (Ha/H/3) and various properties of the galaxies. The meaning of the numbers in the upper right 
is the same as in Fig. 5. The thick curve represents the median value of log (Ha/H/3). The thin curves delimit the zone containing 
80% of the data points (see Sect. 3.1). 
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Fig. 8. Same diagrams as in Fig. 6 except that the ordinate is C(H/?)/([N n] i6584/Ha) u y instead of log (Ua/HB) (see Sect. 3.3) 



hand, the galaxy surface brightness SB is virtually independent 
of eClass (panel b). 
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Fig. 7. The relation between the metallicity indicator [N n] 
/16584/Hq' and the absolute magnitude Mb in the galaxies of 
our sample. Same layout as for the panels in Fig. 6. 

3. Balmer extinction of normal galaxies 

3.1. Definition of the diagnostic plots 

The value of the Balmer extinction is derived from the observed 
ratio of Ha and H/3 line intensities using the relation 



c(m = 



i 



log 



7(Ha)/7(H0) 



/(H0)-/(Hflf) b I°(Ha)/I°(H(3) 



(1) 



where 1° (Ha) /I (H/3) is the intrinsic intensity ratio of these 
two lines and f(H/3) - f(Ha) is equal to 0.335. The value of 
1° (Ha) /I (H/3) is essentially insensitive to the physical con- 
ditions of the gas, being 3.03 at a temperature of 5000 K, 
2.86 at 10000 K, and 2.74 at 20000 K (Osterbrock 
Thus, there is a practically linear relation between C(H/3) 
and log (I(Ha)/I(H/3)) (in the following simply noted log 
(Ha/Hff)). 

In Fig. 6, we plot the value of log (Ha/H/3) as a function 
of various parameters characterizing the galaxies. In all the 
panels, the presentation is the same. The meaning of the num- 
bers reported on the upper right is the same as in Fig. 5. The 
thick curve represents the median value of log (Ha/H/3) and the 
thin curves represent the 10 and 90 percentiles, meaning that 
10% of the data points are respectively above or below these 
curves (the data are binned according to the abscissa in bins 
with equal numbers of points). Such a representation, associ- 
ated with the relevant statistics, is believed to give a fair idea of 
the behaviour of our sample in the different planes. Panel a con- 
cerns the galaxy spectral type eClass, panel b the colour (u — i) 
as obtained from the SDSS data base. Panel c the discontinuity 
at 4000 A, D(4000), defined by the ratio of the flux from 4050 
to 4250 A to the flux from 3750 to 3950 A. Thus, panels a, b, 
and c relate Ha/H/3 to quantities linked to the stellar continuum. 
Panels d, e, and / involve respectively the absolute magnitude 
M(r), the surface brightness SB defined above and the concen- 
tration index CI defined as the ratio of petroR50 to petroR9® 
retrieved from the SDSS data base. Panel g concerns the total 
radius R of the galaxy in kpc, computed from the petroR9® ra- 
dius. Thus, one can consider that panels d, e, f, and g are linked 
to the masses and mass distributions of the galaxies. Panel h 
concerns the inclination parameter b/a derived from the SDSS 



data base. The fact that the SDSS contains such a large num- 
ber of galaxies allows us to infer conclusions from the variation 
of observed properties with galaxy inclination, which was not 
the case in our previous studies (SS99, SS01). Panel i concerns 
UH/3) COTr , the total luminosity of the galaxy in H/3 (in solar lu- 
minosities) after correction for extinction using the value of 
C(H/3). Panel j concerns the equivalent width of Ha. Panels 
i and j are thus related to the importance of H n regions in 
the global galaxy spectra. Finally, panels k and I concern the 
metallicity of the galaxies. Among the various metallicity in- 
dicators, we have chosen two which are complementary. The 
first one, [N n] /16584/Hq' (panel k) has been proposed rather 
recently as a metallicity indicator in giant H n regions (van Zee 
et al. 1 19981 Denicold et al. 2002). Its advantages are that it in- 
creases monotonically with metallicity and is independent of 
reddening, since the two lines it involves have similar wave- 
lengths. The other one, [O in] i5007/[N n] ,16584 (panel was 
first proposed by Alloin et al. ( 1979 ). This indicator also varies 
monotonically with metallicity (it is smaller at large metallic - 
ities) which, in a study of the present kind, makes it safer to 
use than the traditional ([O n] /l3727+[0 m] A5001)/H/3 indi- 
cator first proposed by Pagel et al. ^319 1 and still widely in 
use. The [O m] /15007/[N n] /16584 indicator has one advan- 
tage over the [N n] /16584/Ho' one: its larger dynamical range 
in galaxies. Unfortunately this indicator is affected by redden- 
ing. With the SDSS data, we are able to correct line intensities 
for reddening using the Ha/H/3 ratio, and we use the reddening- 
corrected values in our study, assuming that the intrinsic value 
of Ha/H/3 is 2.9. When using the [O m] ^5007/[N n] ,16584 
or [N n] /16584/Ha line ratios, we rejected the objects where 
the measured equivalent width in any of these lines is smaller 
than 1A, in order to avoid spurious values of the ratios linked 
to large uncertainties in the measurements of line fluxes. 



3.2. Systematics of Balmer extinction in galaxies 

Among the investigated relations between the Balmer extinc- 
tion and various parameters of the galaxies, the most outstand- 
ing correlation is the one between Ha/H/3 and the galaxy spec- 
tral types eClass, for which r$ - -0.665, followed by the one 
with galaxy colour (u - i) for which r$ = 0.630. We find that 
early-type and redder galaxies suffer higher Balmer extinction. 
Of course, without any modelling, one could be tempted to 
say that early-type galaxies are redder not intrinsically, but be- 
cause they suffer larger reddening. Such an extreme explana- 
tion would of course be unreasonable, since there are many ar- 
guments showing that early-type galaxies are indeed composed 
of older stellar populations (e.g. Sommer-Larsen 1996 Fioc 
& Rocca-Volmerange 1997 1. There are, however, indicators of 
the characteristics of the stellar population that do not depend 
on reddening, like for example the discontinuity at 4000 A, 
Z)(4000). This discontinuity is larger for older and for more 
metal-rich stellar populations. This is due to higher Balmer 
opacity for lower stellar effective temperatures dominating the 
spectra in the first hypothesis, and higher metal opacity in the 
second case. It has been argued on empirical grounds (Dressier 
& Shectman 1987 1 that D(4000) is rather sensitive to star for- 



8 



G. Stasiriska et al.: What drives the opacity of galaxies? 



t oo 

en 
O 



9733 
-0.501 
0.0E+00 
-0.691 




-i— 




-22 -20 -16 -12 -20 

M(rJ [mag] M(r) [mag] 

Fig. 9. L(H/?) COIT versus M(r) (left) and D(4000) vesus M(r). Same layout as Fig. 5. 



-16 



mation, but others (e.g. Poggianti & Barbara 1997 1 point out 
that metallicity is important. In any case, the correlation we 
observe with D(4000) (r s - 0.501) indicates that the Balmer 
extinction is higher for galaxies with older and/or more metal- 
rich stellar populations. That the correlation is slightly worse 
than in panels a and b probably reflects the fact that D(4000) is 
not sensitive to reddening. 

We also see a strong correlation with M(r) (rs = -0.45 1) ex- 
pressing that the Balmer extinction is higher for more luminous 
(thus likely more massive) galaxies. Correlation with other pa- 
rameters linked with galaxy stellar mass and concentration are 
present but are less important (rs - -0.298 for SB, rs = -0.182 
for a and r s =0.182 for R). 

Surprisingly, there is no correlation with b/a [rs = -0.091), 
while one would expect the extinction to be larger for more in- 
clined galaxies. One interpretation might be that most of the 
Balmer extinction comes from dusty blobs which are more 
or less uniformly distributed in the galaxy, so that the ex- 
tinction of the integrated flux does not depend on inclination. 
However, such an interpretation is not compatible with the fact 
the near-infrared colours of spiral galaxies depend on inclina- 
tion (Masters et al. 2003). We will come back to this point later. 

There is a weak anticorrelation with EW(Ha) (rs = -0. 194). 
EW(Ha) is a star formation history indicator (Kennicutt 1998) 
and the observed anticorrelation suggests that the Balmer ex- 
tinction is higher in galaxies with a larger proportion of old stel- 
lar populations. On the other hand, the correlation with the total 
HB luminosity (rs = 0.607) suggests that the Balmer extinction 
is larger for galaxies with a larger total amount of young stars. 

Finally, correlations with the metallicity indicators are 
fairly important (panels k and 0: r s = 0.418 for [N n] ^6584/Ha 
and r s = -0.449 for [O m] ^5007/[N n] i6584, suggesting that 
metallicity and Balmer extinction are strongly correlated. 

3.3. Is the metallicity of galaxies the only driver of 
Balmer extinction? 

It is known that many properties of galaxies are correlated in 
the Hubble sequence of galaxies (Roberts & Haynes 1994 1. 
Early-type galaxies tend to be more massive, composed of 
older stellar populations and more metal-rich (see e.g. Zaritsky 
et al. lT994l Kennicutt Il998l Fioc & Rocca-Volmerange[l997). 



Could it be that the main parameter responsible for the trends 
mentioned above between C(HB) and other galaxy parameters 
is simply metallicity? 

Let us assume for simplicity that C(HB) can be written as: 

C(UB) = f(Z)g(X), (2) 

where f(Z) is a function of metallicity and g(X) is a function 
of other parameters. 

We know that [N n] ^6584/Ha and [O in] ^5007/[N n] 
/16584 are metallicity indicators and that they vary monoton- 
ically with metallicity. The problem is that we do not know the 
exact relation between the value of these ratios and the average 
metallicity of the interstellar medium in the galaxy. It is not 
possible to directly use the calibrations proposed for giant H n 
regions (which in themselves are already a matter of debate, 
see Pilyugin 2003 1 because, as mentioned by SS01, spectra of 
galaxies are affected by metallicity gradients and by the con- 
tribution of diffuse ionized regions to the spectrum. This can 
be illustrated as follows. In the sample of spiral galaxies stud- 
ied by Zaritsky et al. ( 1994) the average metallicities vary from 
12 + log (O/H) = 8.3 to 9.3 when the galaxy absolute B mag- 
nitudes Mb range from -17 to -21. The slope of the relation 
between log O/H and Mb in their Fig. 10 is roughly -0.25 (but 
there is substantial dispersion). If we would apply the calibra- 
tion of Van Zee et al. JT9981 ( 12 + log (O/H) = 1 .02 log ([N n] 
^6584/Ha) + 9.35), the corresponding [N n] ^6584/Ha values 
would range between -1.04 and -0.06 and the slope between 
log [N n] /16584/Ha' and Mg would again roughly be -0.25. 
Figure 7 shows the relation between log ([N n] /16584/Ha) and 
Mb in our sample (the value of Mb has been estimated using 
the transformations between the SDSS photometric system and 
the standard UBV system following Smith et al. 2002 1. Most 
of the observed values of log([N n] /16584/Ha) are between 
-0.7 and -0.3 even for the galaxies with the lowest luminosi- 
ties in our sample. This is likely due to the fact that in spectra of 
galaxies both the radial metallicity gradients and the contribu- 
tion of emission from the diffuse ionized medium would tend 
to enhance the [N n] /i6584/Ha ratio. 

What we can do, however, is study the function 
C(H/?)/([N n] /16584/Hq') £I where a is determined empirically 
so as to remove any dependence of this function on [N n] 
^6584/Ha, implying that C(BB)/([N n] ^6584/Ha) fl is inde- 
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pendent of metallicity. We find that a is approximately equal 
to 0.9 in our sample. While this procedure is very schematic 
and relies on assumptions, the fact that we find C(H/3)/([N n] 
/}6584/Ha) 9 to be independent of the other metallicity indi- 
cator, [O m] /15007/[N n] /16584, gives us confidence in our 
approach. (Note that the value of a depends on the definition 
of our sample and should not be used for other samples with- 
out caution). Figure 8 shows the same plots as Fig. 6, with log 
(Ha/HjS) now replaced by C(H0)/([N n] i6584/Ha) 09 . We see 
that C(H/?)/([N n] ^6584/Ha) 9 (which can be assimilated to 
g(X) in Eq. 2) is strongly correlated with the galaxy spectral 
type eClass (r s = -0.436) and colour (u - i) (r s = 0.394) 
and also with D(4000) (r s =0.340) although less tightly than 
C(H/3). The strongest correlation in this figure is with L(H/3) colT 
fa = 0.512). However, we show in Fig. 9, that L(H/3) con is very 
strongly correlated with M(r) (r s = -0.691) so that we conclude 
that an important driver of the Balmer extinction is actually the 
total galaxy luminosity (or mass). On the other hand, D(4000) 
is not correlated with M(r) as seen in Fig. 9, so the observed 
correlation between C(H/3)/([N n] ^6584/Ha) 9 and D(4000) 
implies that the presence of old stars is in itself also a driver of 
the Balmer extinction. 

We also note in Fig. 8 (panel h) that C(H/3)/([N n] 
A6584/Ha) 0S> is linked with galaxy inclination (r s = -0.200). 
One explanation that can account for this as well as for the 
absence of correlation between HafH/3 and b/a is that H n re- 
gions are rather close to the galactic plane of galaxies (which 
is in better agreement with the observed height scale of H n 
regions in nearby galaxies), and that the increased path length 
of the attenuating dust in inclined galaxies is roughly compen- 
sated by abundance gradients: when a galaxy is inclined, the 
innermost regions, that are more metal-rich and therefore most 
likely dust-rich, are more obscured than the outermost H n re- 
gions, which suffer less extinction at optical wavelengths. The 
compensating effect of the metallicity gradients at near-infrared 
wavelengths is far less important due to smaller optical depths 
at these wavelengths, which explains why Masters et al. (2003) 
do find a correlation between near-infrared colours and inclina- 
tion. We have tested by simple toy models that such an expla- 
nation can be viable. 

4. Discussion 

We have thus found (confirming the results of SS99 and SS01) 
that the Balmer extinction of galaxies decreases steadily from 
early- to late-type spirals (whether classified according to their 
spectral type or according to their colour). We have also shown 
that there is a direct dependence of Balmer extinction on metal- 
licity, and also on the age of the stellar populations and on the 
total galaxy luminosity. This last statement is consistent with 
the finding of Wang & Heckman (1996), based on relations 
using far-ultraviolet and far-infrared fluxes, that "the optical 
depth of normal galactic disks increases with galaxy luminos- 
ity". On the other hand, the dependence we find on galaxy type 
seems to contradict the general opinion summarized by Calzetti 
(2001 1 and the recent result by Kauffmann et al. (2003) who 
conclude that "galaxies with the youngest stellar populations 
are the most attenuated by dust". While previous studies were 
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Fig. 10. Histogram of the values of the Balmer extinction at H/3, 
C(HyS) in our sample of galaxies. 

sometimes based on small samples, the result of Kauffmann 
et al. (2003) comes from a model fitting of the continua of 10 5 
galaxies from the SDSS. It should be noted that what they mea- 
sure with this elaborate procedure is the attenuation of the stel- 
lar light. Our work concerns the extinction of nebular light and 
the determination of the Balmer extinction is straightforward. 
The interpretation in terms of global opacity of galaxies is how- 
ever not simple: with Hc/HjS we measure some sort of average 
opacity of the zones dimming the light from H n regions (ex- 
cluding the most opaque ones). Still, it is striking to see how 
well C(HjS) correlates with the galaxy spectral type. In this sec- 
tion, we present a possible phenomenological interpretation of 
our results, and show that our findings can actually be recon- 
ciled with the results of Kauffmann et al. (2003). 

4.1. A phenomenological interpretation of the 

observed Balmer extinction trends in our sample 



Let us write C(H/3) as the following product: 



C(H{3) = n d lcr d 



»d "m "h 



(3) 



where «a is average number of dust particles per unit volume in 
a galaxy, « M is average number of metallic atoms (in any form) 
per unit volume, n H is the average number of hydrogen parti- 
cles per unit volume, n„ is the average number of stars per unit 
volume, / is the geometrical thickness of the region responsible 
for the optical extinction, and cr d is the typical extinction cross 
section of dust grains. 

The common view that the opacity should increase from 
early- to late-types is comforted by such arguments that "galax- 
ies with young stars contain more gas and hence more dust than 
galaxies with old stellar populations" (Kauffmann et al. 2003 1. 
On the other hand, early-type galaxies are more metal-rich and 
on average more massive than late-type galaxies, so it is not 
necessarily surprising that the extinction is actually larger for 
late-type galaxies. In Sect. 3.2 we have given arguments to say 
that the Balmer extinction is also determined by the mean age 
of the stellar population. This could indicate that in late-type 
galaxies «d/«M is larger than in early-type ones. Such a view 
seems indeed to be supported by recent models for the evolu- 
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tion of dust in galaxies (Hirashita 1999 1 which take into ac- 
count the processes of formation and destruction of dust (note 
however that there are presently many assumptions and uncer- 
tainties in such models, see e.g. Dwek 1998 or Edmunds 2001 1. 
In these models, condensation in cool stellar winds from low- 
mass stars is an important source of dust production. 

We note that 75% of the galaxies of our sample have a mea- 
sured C(H/3) between 0.3 and 0.9, as seen from the histogram 
shown in Fig. 10. On the other hand, judging from the metallic- 
ities derived for galaxies of luminosities similar to those of our 
sample (Zaritsky et al. 119941 Chariot et al. I2002i . the metal- 
licity range in our sample is likely higher than just a factor of 
three. This suggests that factors other than metallicity act to re- 
duce the observed range in C(H/3). Obviously, «h/«* is a good 
candidate, since it decreases from late- to early-type galaxies 
(Roberts & Havnes [T994l 

On the other hand nj, which can be assimilated to the stel- 
lar surface density, decreases from early- to late-types. 

The last factor in Eq. (3) is <x d , and one might expect some 
systematic effects if the grain size distribution depends on the 
processes for grain growth or destruction that could have dif- 
ferent relative importances in galaxies of different types. 

In conclusion, dust extinction in galaxies involves many 
factors, and our finding that C(H/3) increases from late- to 
early-types can easily be accounted for, at least qualitatively, 
within our present-day understanding of galaxies and their con- 
stituents. 

4.2. The extinction of stellar light versus the extinction 
of nebular light 

As emphasized above, what we measure with C(H/3) is actually 
the reddening of the nebular light (called Ci in SS01). On the 
other hand, what is determined by Kauffmann et al. (2003) is 
the reddening of the stellar light (called C c in SS01). We now 
examine the relation between these two quantities in our sam- 
ple of galaxies. 

It had already been noted by SS99 and SS01, on much 
smaller samples of galaxies, that EW(Ha) and EW(Hy6) corre- 
late extremely well (see also Kennicutt 1992). This implies that 



the difference between the extinction of the stellar continuum 
and that of the nebular emission is strongly linked to the colours 
of the galaxies, being larger for redder galaxies, (see Eq. 5 in 
SS01). In Fig. 1 1 we show the values of EW(H/8) as a function 
of EW(Ha) in our sample of galaxies. The correlation is ex- 
tremely strong, with tk = 0.782 and r$ = 0.93 1 . Assuming that 
the relation between these emission line equivalent widths is 
EW(H/?) = Ax EW(Ha) (the same model as adopted in SS01), 
an ordinary least square bisector linear fitting (Isobe et al. 1990) 
gives A = 0.185 + 0.001. Considering only objects for which 
EW(Ha) > 20 A the result is the same, A = 0.185 + 0.001. 
These results may be compared with A = 0.194 + 0.011 and 
A = 0.245 + 0.007, presented in SS99 and SS01, respectively. 
Our value A is consistent with that obtained by SS99 but not 
with that obtained by SS01 . The cause of this discrepancy is not 
clear; it might be due to differences in the spectral resolution, 
because the spectral resolution of the SDSS spectra is almost 
half that of the NFGS spectra (Jansen et al. l2000all2000bl used 
bySSOl. 

From Eq. (5) in SS01 and the value of A, we derive that 



C c = C, -0.81 -2.99 log 



Fg(Ha) 



(4) 



where F"(Ha) and F°(H/3) are the intrinsic (i.e. not affected by 
extinction) stellar fluxes in the continuum adjacent to Ha and 
H/3 respectively. We can use the dust-free spectrophotometric 
models of Barbara & Poggianti (1997) to relate the value of 
F°(Ha)/F°(Hf3) to the values of D(4000) and EW(Ha). Using 
panel j of our Fig. 6 we can relate EW(Hff) to log (Ha/Hfi) 
(by taking the median value of Ha/Hfi for a given EW(Ha)) 
and therefore to C/. With the help of Eq. (4) we thus find that 
C c goes from about 0.02 for D(4000) = 1.7 to about 0.2 for 
Z)(4000) = 1.3. This goes in the same direction as the results of 
Kauffmann et al. (2003) who find that the median value of A,, 
which is approximately equal to our C c , roughly goes from 0.2 
at £>(4000) = 1 .7 to 0.6 at D(4000) = 1.3. That the numbers are 
not exactly the same as the ones found by our analysis is not 
necessarily a worry given the dispersion in the observational 
points (both here and in the work of Kauffman et al. 2003) and 
given that the definition of the derived extinction is not exactly 
the same. 

As for C;, it can be evaluated using Fig. 6c and Eq. 1. We 
find C, 0.8 for D(4000) = 1.7 and Q * 0.3 forD(4000) = 1.3 
We note that C c is smaller than C/ at both extremes of the spec- 
tral type range, which can be interpreted as due to the fact that 
dust is more concentrated (and thus more opaque to radiation) 
in molecular clouds associated with H n regions than in the dif- 
fuse interstellar medium. One may wonder why C c increases 
from early- to late-type spirals while C/ decreases. The answer 
to this may be related to the fact that the stellar light from early- 
type galaxies is dominated by the bulge and to a specific dis- 
tribution of dust resulting from dynamical effects. Advanced 
3D-modelling of the star, dust and gas distribution in galaxies 
would be needed to test any interpretation of our empirical re- 
sult. 
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Fig. 12. Galaxies from our sample detected by IRAS (see Sect. 5). Plots of log (Ha/Hfi) vs. log L(FIR) (panel a) and log L(H/3) 
(corrected for Balmer extinction) vs log L(FIR) (panel b). 



5. The relation between C(H/3) and infrared 
luminosity in our sample 

Another indicator of the presence of dust in galaxies is their 
far-infrared radiation attributed to the emission from dust 
grains heated by stellar radiation (e.g. Inoue 2002). Kewley 
et al. (2002), using 81 galaxies with relevant data from the 
NFGS sample, noted the existence of a very strong correlation 
between the far-infrared luminosity L(FIR) and the Ho- lumi- 
nosity (corrected for reddening using the observed HafH/5 ra- 
tio). They find a Spearman rank correlation coefficient of r$ = 

0. 98. They do not plot the relation between HafHfi and L(FIR), 
but this plot is shown by Wang & Heckman ( 1996 1 for about 
30 galaxies from the Kennicutt ( 1992 > atlas with relevant data, 
and they find a correlation to within better than a 0. 1 % level 
of confidence (the correlation coefficient from a standard linear 
regression is 0.63). 

It was tempting to produce such diagrams for our sample of 
galaxies using the IRAS data base. We have searched the IRAS 
Faint Source Catalogue (Moshir et al. 119891 and the IRAS 
Point Source Catalogue (IPAC 1986 1 for far-infrared counter- 
parts of galaxies from our sample. Since the IRAS beam size at 
60 fim is V.5 and the IRAS uncertainty in position is about 30", 
we used a conservative circular window of 1' radius around 
the SDSS positions of our sample galaxies to search for far- 
infrared sources. Note that we also constituted a subsample 
of IRAS galaxies using a search window of 30" instead of 1'; 
this subsample is twice smaller but the conclusions drawn be- 
low are not altered. Since an estimate of the total infrared lu- 
minosity requires IRAS fluxes at both 60 fim and 100/im (fa 
and /ioo, respectively), we kept only those galaxies with avail- 
able fluxes in both bands and IRAS quality flags greater than 

1. About 80 galaxies were common to both the IRAS Point 
Source Catalogue and the IRAS Faint Source Catalogue. For 
those objects, we retained the flux values given by the Faint 
Source Catalogue, which have greater quality. 

This search resulted in 189 galaxies. We have computed 
L(FIR) using the formula from Thuan & Sauvage ( 119921 : 

L(FIR) = 3.95 x 10 5 (2.58/ 60 + fao)D 2 h (5) 

where L(FIR) is in solar units, fa and /ioo are in Jy and is 
the cosmological luminosity distance of the galaxy in Mpc. In 



most cases, L(FIR) > 2 x 10 10 L . Note that this lower limit 
actually corresponds to the high luminosity tail of the lumi- 
nosities of spiral galaxies detected by IRAS (Thuan & Sauvage 
1992). This stems from the fact that our samples contains only 
galaxies at redshifts larger than 0.05, and at such distances only 
luminous infrared galaxies can be detected. 

Figure 12 plots the values of log(Ha/H/3) (panel a) and 
log L(HyS) C01T (panel b) versus logL(FIR). We do not see any 
correlation between log (Ha/H/?) and L(FIR) (the probability 
associated to the Kendal statistics is as large as 0.2). A corre- 
lation between L(Hy8) con and L(FIR) is present in our sample, 
but much weaker than found by Kewley et al. (2002). It thus 
turns out that correlations between optical and infrared prop- 
erties become much less conspicuous or even disappear when 
restricting the sample to luminous infrared galaxies 1 . Note that 
in the samples of Wang & Heckman ( 1996 1 and Kewley et al. 
(120021 galaxies with L(FIR) > 2 x 10 10 L Q are only a handful 
and do not show a clear correlation either. 

It is also interesting to see how our luminous infrared galax- 
ies compare with our entire sample of galaxies. For this, we 
show in Fig. 13 the same diagrams as in Fig. 6, with the lu- 
minous infrared galaxies now indicated by large circles. We 
see that in all the diagrams the galaxies detected by IRAS are 
found in the entire domain covered by the SDSS galaxies ex- 
cept the regions corresponding to the least luminous galaxies 
(which cannot be detected by IRAS at the redshits of our SDSS 
sample). The distributions of the luminous infra-red galaxies is 
very similar to that of our entire set of SDSS galaxies, except 
that is is slightly skewed towards the highest values of Ha/H/?. 

The comparison of our SDSS sample with IRAS data raises 
two questions. Why were only 2% of our sample of SDSS 
galaxies detected by IRAS? Why are those IRAS detected 
galaxies found at any value of Hc/H/?? Concerning the first 
question, we have already noted that, because of the limit im- 
posed on the redshift in our SDSS sample, only luminous in- 
frared galaxies can be identified as IRAS sources. An IRAS 
source is characterized by the presence of warm dust (50- 
200 K), which lies relatively close to the stars. Such a situa- 

1 We use the term "luminous infrared galaxies" for convenience, 
although it is traditionally reserved for galaxies with L(FIR) > 10" 
L (Sanders & Mirabel 1995J. 
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Fig. 13. Same diagrams as in Fig. 6, with the galaxies detected by IRAS represented with large circles. 



tion is encountered in zones where star formation has occured 
very recently, and which have not yet been swept out by stel- 
lar winds. In a galaxy, the H n regions that emit the observed 
H B aimer lines have ages between 1 and lOMyr, but the H n 
regions that will emit most efficiently in the IRAS bands are 



younger than that. This explains, at least qualitatively, why only 
a small proportion of the galaxies from our SDSS sample have 
been detected by IRAS. These are the ones that happen to have 
a larger number of zones of very recent star formation. Such 
galaxies are more likely found among more massive galaxies at 
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a given spectral type, but may be found at any galaxy spectral 
type. Since the IRAS emission probes these extremely young 
star forming regions, there is no reason why it should show any 
correlation with the global optical properties of the galaxies, 
which answers the second question. 

6. Summary, implications and prospects 

We have used the observations from the First Data Release of 
the SDSS to examine a sample of normal galaxies (as opposed 
to galaxies with an active nucleus) and to investigate how the 
Balmer extinction C(HB) (i.e. the extinction at the wavelength 
of H/3 derived from the HafH/3 emission line ratios) relates with 
other global properties of the galaxies. Our selection criteria to 
build up the sample resulted in a data set of 9840 galaxies with 
adequate data. All these galaxies are at redhifts larger than 0.05 
to avoid strong aperture effects. 

Our main findings are the following: 

1. C(HB) is linked with the galaxy spectral type and colour, 
decreasing from early- to late-type spirals. 

2. C(HB) increases with increasing metallicity 

3. C(HB) is, probably, also affected by the age of the stellar 
population, being larger in the case of older stellar popula- 
tions. 

4. C(HB) depends on galaxy masses. 

5. The extinction of the stellar light is correlated with both 
the extinction of the nebular light and the intrinsic galaxy 
colours, resulting in a trend with galaxy colour that may be 
opposite to the trend of C(HB). 

The present work thus confirms the conclusions of our pre- 
vious studies (Sodre & Stasinska 1999 and Stasinska & Sodre 
2001 1, which were based on much smaller samples and used 
data with lower spectral resolution (Kewley et al. 2002 using 
the same sample as SS01, also found that early-types in that 
sample are more heavily reddened than late-types). 

Compared to our previous studies, the large number of 
galaxies in the SDSS sample allows us to investigate issues re- 
lated to the inclination of galaxies. 

The fact that C(HB) correlates so well with other properties 
of galaxies is remarkable, given that the extinction, especially 
in late-types, is known to be not uniform across the face of 
galaxies (e.g. Beckman et al. 1996 1. 

We have cross-correlated our sample of SDSS galaxies with 
the IRAS data base in order to investigate any relationship be- 
tween C(HB) and total infrared luminosity of the galaxies. Due 
to the lower redshift limit imposed to our sample and to the de- 
tection limit of IRAS, such a procedure selected only luminous 
infrared galaxies. We found that correlations that were shown 
by other authors to exist between optical and infrared proper- 
ties of galaxies disappear when restricting to luminous infrared 
galaxies. We also found that the optical properties of the lumi- 
nous infrared galaxies in our SDSS sample are very similar to 
those of our entire sample of SDSS galaxies. 

We have proposed a phenomenological interpretation of 
our findings. We suggest that the main driver of the Balmer 
extinction of galaxies is their mass, combined with their metal- 
licity and presence of old stellar populations. The infrared 



luminosity of the galaxies as determined by IRAS, which 
is attributed to radiation from hot stars reprocessed by dust 
grains, samples the regions with the most recent episodes of 
star formation, and is not connected with the Balmer extinc- 
tion. Obviously, detailed modelling of the spectral light from 
galaxies taking into account the effects of dust and using a 
complete code such as GRASIL (Silva et al. 1998, see also 
http://web.pd.astro.it/granato/grasil/grasil.htmli is needed for a 
deeper understanding of the empirical relations we have found. 
This is not an easy task, however, since as noted by Witt at al. 
(1992) and Witt & Gordon (2000), equal amounts of dust in 
different configurations may produce very different reddening 
and attenuation effects. In any case, future models of the inte- 
grated light from galaxies including the effects of dust should 
also aim at reproducing the correlations we have shown. 

An important outcome of our study is to open the way for 
an improved correction for extinction in the determination of 
such parameters as the global star formation rate in galaxies or 
their total stellar masses. For normal galaxies, the global star 
formation rate can be obtained from the total Ha luminosity 
corrected for extinction using the extinction derived from the 
HafH/3 emission line ratio (keeping in mind the reservations 
expressed e.g. by Hirashita et al. 2003). If observations do not 
allow one to determine the Balmer extinction, one can make use 
of e.g. the observed galaxy colour or the D(4000) parameter 
to obtain an estimate of statistical value since we have shown 
that all these quantities are correlated. On the contrary, the total 
stellar mass can be estimated from the observed stellar fluxes of 
the galaxy after correcting for stellar extinction. This should be 
done with a proper model fitting of the observed continuum as 
in Kauffmann et al. (2003). However, the strong correlation that 
we have found empirically between stellar extinction, Balmer 
extinction and galaxy colours can provide a basis for a statis- 
tical method to determine the total masses of galaxies. These 
aspects will be developed in future work and should be impor- 
tant especially for the study of galaxies at intermediate and high 
redshifts. 
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